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It has long been known that particles with short-range repulsive interactions in spatial dimension 
d = 1 form universal quantum liquids in the low density limit: all properties can be related to 
those of the spinless free Fermi gas. Previous renormalization group (RG) analyses demonstrated 
that this universality is described by an RG fixed point, infrared stable for d < 2, of the zero 
density gas. We show that for d > 2 the same fixed point describes the universal properties of 
particles with short-range attractive interactions near a Feshbach resonance; the fixed point is now 
infrared unstable, and the relevant perturbation is the detuning of the resonance. Some exponents 
are determined exactly, and the same expansion in powers of (d — 2) applies for scaling functions 
for d < 2 and d > 2. A separate exact RG analysis of a field theory of the particles coupled to 
'molecules' finds an alternative description of the same fixed point, with identical exponents; this 
approach yields a (4 — d) expansion which agrees with the recent results of Nishida and Son (Phys. 
Rev. Lett. 97, 050403 (2006)). The existence of the RG fixed point implies a universal phase 
diagram as a function of density, temperature, population imbalance, and detuning; in particular, 
this applies to the BEC-BCS crossover of fermions with s-wave pairing. Our results open the way 
towards computation of these universal properties using the standard field-theoretic techniques of 
critical phenomena, along with a systematic analysis of corrections to universality. We also propose 
a 1/N expansion (based upon models with Sp(2iV) symmetry) of the fixed point and its vicinity, 
and use it to obtain results for the phase diagram. 



I. INTRODUCTION 



The study of ultracold atomic gases has drawn renewed 
attention to interacting quantum gases in regimes where 
their properties are independent of all microscopic details 
of the interaction between the atoms. This happens when 
the scattering cross-section between two atoms reaches 
its unitarity limit. 

The simplest example of this phenomenon is the di- 
lute Bose gas with short-range repulsive interactions in 
spatial dimension d = 1 where, in fact, the properties 
are generically universal. It is a well-known property 
of quantum-mechanical scattering in one dimension that 
for any generic short-range potential, the S matrix for 
two particle scattering reaches its unitarity limit value 
of —1 in the limit of zero momentum transfer. However, 
the strongly interacting Bose gas with S = — 1 can be 
reinterpreted^ as a free Fermi gas (the 'Tonks' gas), al- 
lowing easy computation of at least the thermodynamic 
properties. Correlation functions are much harder to 
compute, especially at finite temperature and unequal 
times, but much theoretical effort has been expended in 
this direction^^ Experimental studies of Tonks gas 
behavior have also appeared^ 

No generic unitarity limit scattering is obtained in 
d = 3. In this case, we usually have the limit 5 = 1 for 
low momentum, implying that sufficiently slowly mov- 
ing particles do not scatter. However, upon fine-tuning 
the strength of an attractive interaction, it is possible 
to obtain unitarity limit scattering across a Feshbach 
resonance Near a Feshbach resonance, the s-wave 



scattering amplitude, /o, takes the form 10 



where k is the momentum transfer, S = 1 + 2ikfo, and 
v is the 'detuning' across the resonance (it is related to 
the scattering length, a, by v = — 1/a). The resonance 
occurs at v — 0, and notice then that S = — 1, the uni- 
tarity limit. We will be interested here in broad Feshbach 
resonances, where Eq. ([T]) holds over all the momenta of 
interest in a quantum gas. For a gas with density n, a 
resonance is broad if the 0(k 2 ) corrections to the de- 
nominator of Eq. (JTJ) are negligible at the characteristic 
momentum n 1 / 3, (we use h = 1 throughout). Notice then 
that the width condition on the Feshbach resonance, is 
actually a dilutcncss condition on the quantum gas, and 
any Feshbach resonance is 'broad' for a sufficiently dilute 
gas. 

This paper demonstrates that a unified understanding 
of the distinct universal properties of the quantum gases 
in d — 1 and d — 3 can be obtained in a field-theoretic 
renormalization group (RG) analysis. The RG approach 
was used previousl y 3 ' 5 ' 11 ! 12 for the repulsive Bose gas in 
d < 2, and here we demonstrate that a direct generaliza- 
tion applies to a wide class of quantum liquids in general 
d. We will use the structure of an RG fixed point to argue 
that there is a universal phase diagram as a function of 
relevant perturbations away from the fixed point: these 
perturbations are the temperature, the detuning away 
from the Feshbach resonance, and the density imbalance 
for a two-component quantum gas. Our claims of uni- 
versality are stronger than previous ones which applied 
only at the resonance i^i 9 . In particular, questions on the 
nature of the phases, including the existence of exotic 
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gapless superfluid states, have a unique and universal an- 
swer for a broad Feshbach resonance. We will also obtain 
explicit results for this phase diagram, and for associated 
universal scaling functions, for the two-component Fermi 
gas in d = 3 with s-wave pairing using a 1/N expansion 
of a model with Sp(2iV) symmetry. In particular, we ob- 
tain results on the equation of state for the normal state 
with unbalanced densities and two Fermi surfaces. 

The form of our universal phase diagram has similari- 
ties to work by Sheehy and Radzihovsk y 13 i 14 (and related 
work by others^ 5 -) where they departed from the narrow 
resonance limit in an expansion in the resonance width. 
We discover the same phases and transitions between 
them, except the possible non-uniform phases which are 
not a focus of this paper. The main advance in our work 
is a systematic treatment of the broad resonance limit, 
which is experimentally relevant, and in which univer- 
sality applies. Universality has also been explored using 
a different RG analysis in Ref. they do not give a 
"critical phenomena" interpretation of the results, which 
we claim below yield the most direct interpretation and 
computation tools for identifying universal behavior. 

The previous study 3 -! 5 - ' 11 ' 12 of the dilute Bose gas be- 
gan with an RG analysis of the field theory of the zero 
density gas^ at chemical potential p, = 0, considered as 
a quantum critical point between the vacuum state for 
fi < 0, and the finite density ground state for /i > 0. In- 
teractions between a pair of particles created out of the 
vacuum are characterized by a coupling u (defined more 
precisely in the body of the paper) , and this coupling was 
shown to obey the exact field-theoretic RG equation 



d<2 



du . u 2 
— = (2 - d)u 



(2) 



under a rescaling of length scales by a factor of e e , 
and a dynamic exponent z — 2. As we will review 
later, this RG equation applies equally to two-body in- 
teractions between particles with arbitrary statistics and 
masses. As shown in Fig. [TJ for d < 2, this RG equa- 
tion has an infrared stable fixed point at u* > 0. It was 
argue d 5 ' 11 ' 12 that this fixed point interaction which de- 
scribes unitarity scattering. The chemical potential, /i, is 
a relevant perturbation at this fixed point, and its scal- 
ing dimension, and RG flow can also be determined ex- 
actly. A standard 1 ^ field-theoretic analysis then allowed 
computation 5 , of all universal properties of the \l ^ 
Bose gas in d < 2, in an expansion which can be carried 
out to all orders in u* = 2(2 — d). The degree of difficulty 
of this computation is the same for thermodynamic and 
dynamic correlations. Practically speaking, this expan- 
sion is then useful for those physical properties which are 
difficult to access by the fermionization method in d = 1. 
In some cases, an exact computation is possible both in 
the fermionization approach, and also by resumming all 
orders in the (2 — d) expansion, and exact agreement was 
obtained. 3 A review of this approach to Bose and Fermi 
quantum liquids in d < 2 appears in Section HT1 

Let us now examine the RG equation ^ for d > 2. 
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FIG. 1: The exact RG flow of Eq. ©, as discussed in Ref.0. 
Here u is a measure of the short-range two-body interaction 
between the particles, and the RG applies in the limit of low 
density of either Bose, Fermi, or Bose-Fermi quantum liquids, 
(a) For d < 2, the infrared stable fixed point at u — u* > 
describes quantum liquids of either bosons or fermions with 
repulsive interactions which are generically universal in the 
low density limit. In d = 1 this fixed point is described by 
the spinless free Fermi gas ('Tonks' gas), for all statistics and 
spin of the constituent particles, (b) For d > 2, the infrared 
unstable fixed point at u = u* < describes the Feshbach 
resonance which obtains for the case of attractive interactions. 
The relevant perturbation (u — u*) corresponds to the the 
detuning from the resonant interaction. 



The flows are shown in Fig [T] Now there is an infrared 
unstable fixed point atati* < 0. We show here that this 
fixed point describes the Feshbach resonance of particles 
with attractive interactions. The unitarity of the scat- 
tering at this fixed point has also been pointed out inde- 
pendently in Ref. The relevant perturbation (u — u*) 
is proportional to the detuning, v, of the resonance. By 
an elementary computation of the eigenvalue of linear 
perturbations from the u = u* fixed point of Eq. ((2|) we 
therefore determine the exact scaling dimension 



dim[i/] — d — 2 



(3) 



The universal properties of the quantum liquid with v =^ 
and fi ^ can now be determined by essentially the 
same renormalized theory earlier used for d < 2, and now 
becomes an expansion in (d — 2). Indeed, as we will see 
in Section [TX|, the earlier— (2 — d) expansions for scaling 
functions for d < 2, apply unchanged at resonance for 
d > 2. A shortcoming of this (d — 2) expansion for the 
attractive Fermi gas is that the pairing amplitude is of 
order exp(l/u*) ~ exp(— l/(d — 2)), and so all effects 
associated with superconductivity are non-perturbative. 

Regardless of the ability to determine scaling functions 
in the (d — 2) expansion, the flow equation (J2]), and its 
associated exponents are exact, and allow us to deduce 
some exact scaling forms associated with the RG fixed 
point in d > 2. For the two-component Fermi gas, we 
have already noted two relevant perturbations at the 
fixed point, the detuning v and the chemical potential 
/x. For general d, we define the overall scale of v by the 
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FIG. 2: Universal phase diagram at zero temperature (T = 0) 
and balanced densities (h = 0) for the two-component Fermi 
gas in d = 3. The vacuum state (shown hatched) has no par- 
ticles, and is present for \x < and v > 0, or for v < and 
H < —v 2 /(2m), where m is the mass of a fermion. The posi- 
tion of the v < phase boundary is determined by the energy 
of the two-fermion bound state. The density of particles van- 
ishes continuously at the second order quantum phase transi- 
tion boundary of the superfluid phase, which is indicated by 
the thin continuous line. The quantum multicritical point at 
fj, = v = h = T = (denoted by the filled circle) is the RG 
fixed point which is the basis for the analysis in this paper. 

requirement that the s-wave scattering amplitude is pro- 
portional to \/(y— (— fc 2 )( d ~ 2 )/ 2 ). The scaling dimension 
of /i is exactly 2. There is also a third relevant pertur- 
bation at the fixed point: a 'field' h conjugate to the 
density difference of the two species of fermions, which 
also has scaling dimension 2. From all this information, 
and the dynamic exponent z = 2, we can deduce one of 
our central results, a scaling form for the grand canoni- 
cal free energy density F (from the usual thermodynamic 
identity, F — — P, the pressure) at a temperature T: 

r-W^r^.^-J^), (4, 

where m is the mass of one of the particles, and Fj? is 
a universal scaling function whose three arguments ex- 
tend over all real values (if the mass of the particles are 
distinct, the scaling function will also depend upon the 
ratio of the masses). For the two-component Fermi gas 
this universal scaling function contains all the informa- 
tion on the BEC-BCS crossoveri 20 i 21 i 22 i 23 The existence 
of this scaling form also implies that there is a universal 
phase diagram as a function of fj,, v, h, and T in the low 
density limit (or equivalently, for a broad resonance) . In- 
formation on all the phases and phase transitions in this 
phase diagram is also contained in Eq. ([4|. At h = and 
T = 0, the form of this universal phase diagram can be 
deduced exactly, and is shown in Fig.[5]for d = 3. There is 
a second-order quantum phase transition line separating 



the zero density vacuum from the low density superfluid. 
The RG fixed point is a multicritical point on this line, 
and the whole phase diagram is constructed as a theory 
of relevant perturbations from this point. Note that there 
is no phase transition at finite density within the super- 
fluid, only a smooth crossover between the BCS and BEC 
limits. However, there is a quantum critical point at zero 
detuning (at "resonance" ) on the phase boundary, and 
we use this critical point as a point of departure to de- 
scribe the entire phase diagram. We will present further 
results on this phase diagram in the body of the paper. 

Given the limitations of the (d — 2) expansion for 
obtaining useful approximations for the scaling func- 
tion, the remainder of the paper presents alternative ap- 
proaches to analyzing the fixed point described above. In 
Section IIII1 we examine another field theory commonly 
used to describe a Feshbach resonance: a theory of mixed 
'atoms' and 'molecules', with the 'molecule' correspond- 
ing to the bound state in the vicinity of the Feshbach res- 
onance (the so-called "two-channel" approach 14,24 ). As 
argued above, it is useful to perform an RG analysis in 
the limit of zero density, when the physics reduces to 
the quantum mechanics of few-particle scattering. The 
RG flow equations can now also be determined exactly, 
and are presented in Section IHH They contain a fixed 
point, which, at first glance, appears distinct from that 
of Eq. However, a careful analysis of its universal 
properties, such as scaling dimensions of operators, shows 
that they are identical to the results quoted above. The 
finite density field theory associated with this fixed point 
becomes weakly coupled when (4 — d) is small. Conse- 
quently, a standard renormalized perturbative expansion 
now allows us to obtain universal properties in powers 
of (4 — d). The universal results of the (4 — d) expan- 
sion are equal to those obtained recently in the important 
work of Nishida and Son^, who motivated the expansion 
from different considerations and use a different compu- 
tational scheme. The special role of d = 2 and d — 4 was 
also discussed by Nussinov and Nussinov^ Our method 
streamlines the expansion (4 — d) onto the standard ap- 
proach used extensively for critical phenomena, and also 
opens the way to a systematic analysis of corrections to 
scaling by a consideration of irrelevant perturbations. 

Further considerations on the structure of the (4 — d) 
expansion shows that it also has some limitations in its 
range of validity. For the two component Fermi gas, 
the naive (4 — d) expansion is restricted to being well 
within the superfluid state, or well within the normal 
state with one Fermi surface. A subtle redefinition of 
parameter scale o 25 ! 27 is required to obtain an expansion 
for the 2 Fermi surface normal state. The latter is one 
of the most interesting new regimes discovered in recent 
experiments^^, and so an unrestricted analysis here 
would be useful. 

With this motivation, in Section IIVI we introduce a 
third approach towards analyzing the finite density field 
theory at the RG fixed point of Eq. [5] This is based upon 
the 1/N expansion, where N is the number of particle 
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species. This expansion holds uniformly for general fi, v, 
h, and T. For the two-component Fermi gas, the critical 
fixed point has Sp(2iV) symmetry, and the expansion is 
based upon that considered earlier in Ref. well away 
from the Feshbach resonance. 

We conclude this introduction by briefly discussing 
quantum liquids with attractive interactions in d < 2; 
these are not directly described by the fixed points dis- 
cussed so far. In such systems, the 'atomic' constituents 
will bind to form 'molecules' and the resulting molecules 
could have repulsive interactions, leading eventually to 
the quantum liquid of the type discussed above for d < 2. 
The crossovers associated with this behavior will be de- 
scribed by an RG analysis which is a combination of 
that presented in Sections [TT] and IIIII — we need to explic- 
itly include the interactions between the atoms and the 
molecules (relevant for d < 2), along with 3-point term 
allowing for a molecule to decay into two atoms (relevant 
for d < 4). We will not present the details of such an 
RG analysis here, but note that these considerations are 
entirely consistent with a recent analysis by Gurarie^L in 
d = 1 using Bethe ansatz methods. 



renormalized dimensionless coupling u by 



K 2 ~ d Z 4 

2mSd 



(8) 



where k is a renormalization momentum scale, Sd — 
2/(r(d/2)(4?r) d / 2 ) is the standard phase space factor, and 
Z4 is the renormalization constant. The ladder diagrams 
lead to the exact result 
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2-d 



(9) 



in the minimal subtraction scheme. From this, we obtain 
the exact RG (3 function for u presented earlier in Eq. |J2]), 
upon rescaling k by a factor of e £ . 

Now we include relevant operators away from this fixed 
point (apart from (it — it*) for d > 2). For the single- 
component Bose gas, there is only one such term, the 
chemical potential 



bp 



drd d x< —fAiJj^ijj 



(10) 



II. RG FOR A FIELD THEORY OF ATOMS 

We begin by considering the theory of the dilute Bose 
gas, from the perspective of a zero density quantum crit- 
ical point, as discussed in Ref. HI We consider a zero 
density theory at the chemical potential /i = 0. This can 
be viewed as a quantum critical point between the finite 
density phase with [i > 0, and the empty vacuum ground 
state for /i < 0. The critical theory of this quantum 
critical point of a boson, ip, is 



S. 



drd d xl V> f 



d_ 

7h 



2m 



+ y ( 5 ) 



It easy to see that an insertion of S p acquires no operator 
renormalization at the zero density critical point, and so 
from Eq. ([7]) we conclude that 



dim[/i] = 2. 



(11) 



For the two-component Fermi gas, we have an additional 
relevant perturbation corresponding to the difference in 
chemical potential of the two species: 



Sf P = 



dTd d x< 



- ("WW +'02V'2) 



(12) 



The RG analysis of the zero density theory does not 
depend upon the statistics of the particles, so we also 
consider in parallel a theory of two-component fermions, 
■0O-, f = 1,2. The critical theory is 



5/ = / drd d 



% + u ip\i>iip2ipi M6) 



The single particle self-energy is identically zero at 
T = because the vacuum state contains no particles. 
Consequently, we can perform an RG transformation 
with the exact dynamic critical exponent z = 2, and 
obtain the exact scaling dimension 



dim[^i] = dimmer] = d/2. 



(7) 



Despite the trivial nature of the vacuum at T = 0, 
there is a non-trivial RG transformation of the two body 
interaction u Q . This comes entirely from the ladder di- 
agram of repeated two-particle scattering. Defining a 



and again we have 



dim[/i] 



(13) 



Now let us examine the universal properties of the fi- 
nite density /i ^ theory in the vicinity of the fixed point 
of Eq. ([2]). We will consider the cases d < 2 and d > 2 
separately in the following subsections. 



A. d < 2 

The fixed point of Eq. ([2]) has repulsive interactions 
with u* = 0. 

The properties of the Bose ga s at this fixed point were 
already considered in Refs. [Blflll. As an example, we note 
the grand canonical free energy density at T = 0, i.e. the 
pressure P. From Eq. (|U), we deduce that this obeys 



P = C h d ii(2mn) d/2 



(14) 
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where C\ is a universal number. A (2 — d) expansion for 
C\ was presented in Ref. la, and the leading result is^ 



C\ - S d 



1 



4(2 - d) 
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+ 0(2 - d) 



(15) 



The corresponding universal properties of the two- 
component Fermi gas at the repulsive u — u* are also 
easy to work out. For simplicity, we present results only 
at T — and h = 0. Then, the result {S} still applies, 
but with a different universal constant C d . This is easily 
computed from the perturbation theory for the pressure 
of the two-component Fermi gas to the first order in uq: 



P = 
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Jo 


(2n) d \2m 




- .y/Smfi d d k - 


+uo 


Jo w 



+ O(u 2 ) (16) 



The universal result is obtained by evaluating 
u = u* consistently to order (2 — d). 



C{ = S d 



(2-d) 
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(17) 



Of course, in both the Bose and Fermi cases considered 
in this subsection, we can obtain exact results in d = 1. 
The thermodynamics of both models are equivalent to 
the free, spinless, Fermi gas. Note that even the thermo- 
dynamics of the dilute two-component Fermi gas is given 
by the spinless Fermi gas. From this, we easily obtain 



C\ = C[ 



f 



2 

3^ 



(18) 



of Eq. ([2]). In addition to the 'elementary' particles, ipi,2, 
we will also allow for a 'composite' molecule, $. For gen- 
erality, we will now allow the masses of the particles to 
be possibly distinct, m\ and rri2- By Galilean invariance, 
the mass of the 'composite' particle is mi 4- rri2- The 
statistics of the two ip a can be arbitrary — they can be be 
Fermi-Fermi, Bose-Fermi^S, or Bose-Bose 33 . The 'com- 
posite' particle $ is then, naturally, a boson, fermion, or 
a boson respectively. 

Our analysis will find a 'critical' fixed point describ- 
ing the unitarity limit theory. All exponents and scaling 
dimensions associated with the fixed point will turn out 
to be identical to those obtained for a seemingly distinct 
fixed point in Section [III suggesting that the two theories 
are in fact identical. 

As in Section HH it is useful to begin with an RG analy- 
sis of the zero density 'critical' theory, which is presented 
in Section IIII Al The perturbations away from the criti- 
cal point are considered in Section lill Bl 



A. Zero density critical theory 



The action of the critical theory is now 

d V 2 



S c — J dTd d x^>pl 



dr 2m n 
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v <9t 2(toi + m.2) 
9o ($^M>2 + V>M$) + • • • I ( 1Q ) 



B. d > 2 

The fixed point of Eq. © is now at u* = -2(d-2) < 0. 

The Bose gas with such attractive interactions is un- 
stable, and further interactions are needed to stabilize 
the theory: we will therefore not consider it further here. 

The Fermi gas is stable and has universal properties 
that can be computed for small (d — 2). Indeed, the 
expansion described in Section III Al applies unchanged 
also for d > 2. In particular, the pressure is still given by 
Eq. I|17p. and this agrees with a recent result 27 obtained 
while this paper was in preparation. Of course, as noted 
in Section [H this expansion knows nothing about pairing 
between the fermions, as such effects are^ 7 - exponentially 
small in 1/(2- d). 



III. RG FOR A FIELD THEORY OF ATOMS 
AND MOLECULES 

This section will use the popular "two-channel" formu- 
lation of the Feshbach resonanc e) 14 ' 24 to obtain an alter- 
native field theoretic description of universal fixed point 



Here ro c is a cutoff-dependent bare "mass" term which 
has to be tuned to place the theory at its critical point i.e. 
describe the unitarity limit scattering. In dimensional 
rcgularization, r$ c = 0. A tree- level scaling analysis with 
z = 2 shows that 

dimh/v] = d/2 ; dim[$] = d/2 ; dim[ 5o ] = (4- d)/2 

(20) 

All other interactions are irrelevant for d > 2. The scal- 
ing dimension of go naturally suggests that we perform a 
renormalization group analysis as an expansion in (4— d). 

It turns out that it is possible to perform an RG analy- 
sis exactly, to all orders in (4 — d). Most of the Feynman 
diagrams vanish because of the zero density of the parti- 
cles. The only non-zero renormalization is a wavefunction 
renormalization of the composite particle The wave- 
function renormalization is computed from the $ Green's 
function 

1 

G$(fc, u>) 



-ilu + k 2 / (2(mi + to 2 )) + roc — S$(fc, lu) 

(21) 

For the above theory, at T = 0, the exact expression for 
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the self energy is 
E$(fc,w) = gl 



dQ 
2^ 



d d p 
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(2w) d (-iQ+p 2 /(2m 2 )) 
1 



H(w-f!) + (fc-p) 2 /(2mi)) 



2 r(i-d/2) 

( 47r )d/2 
X — ILO + 



2m\m 2 
(mi + m 2 ) 

fc 2 



d/2 



2(mi + m 2 ) 



d/2-l 



(22) 



In the last equation, we have dropped a constant cut- 
off dependent term which cancels against r$ c (additional 
subtractions are necessary for d > 4). We now define a 
renormalized coupling g by 



5o = K (4 - d)/2 



2mim 2 



(mi + m 2 ) 



-d/4 



1 



V z$Sd 
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(23) 



where Z$ allows for a non-trivial wavefunction renormal- 
ization of the $ field. It is easy to check that additional 
self-energy and vertex corrections vanish in the zero den- 
sity theory at T = 0, and so there is no wavefunction 
renormalization for the tpa-, nor an independent coupling 
constant renormalization of go. In terms of g, we can 
expand the self energy to leading order in (4 — d), and 
obtain 



£$(fc, oj) 



Z^{A-d) 



2 (mi + m 2 ) 



(24) 



Cancelling poles in the minimal subtraction scheme^, we 
obtain the exact value of Z$: 



the Feshbach resonance. The attractive flow towards the 
g = g* fixed point above corresponds to the 0(k 2 ) cor- 
rections to the scattering amplitude we discussed below 
Eq. |T]): the system can be assumed to be exactly at 
g = g* for a broad Feshbach resonance, while the irrel- 
evant flows towards g — g* have to be included for a 
narrow resonance. As stated earlier, we restrict our com- 
putations here to the broad resonance. 

The critical exponents are to be evaluated at the fixed 
point of Eq. (|26|) . From this, we obtain the exact results 



V = V{9*) = (4 - d) ; dim[$] = (d + 77)/2 = 2. (28) 

This (along with Galilean invariance) implies that the $ 
Green's function behaves like GJ 1 ~ [— iu; + k 2 /(2(mi + 
«i 2 ))] 1 ~' ) / 2 , which is, of course, the result already ob- 
tained in Eq. (gg). 

At this point, the above RG fixed point appears to 
bear little direct relation to the fixed point discussed in 
Section [TTJ In the following subsection we show that the 
spectrum of eigenoperators of the two fixed points are 
identical, indicating that they describe the same univer- 
sal physics. 



B. Perturbations to finite density 

Having obtained the critical theory S c , we can now look 
at all possible perturbations, and classify them accord- 
ing to their renormalization group eigenvalues. It turns 
out that there are 3 relevant perturbations, and all other 
perturbations are irrelevant. Anticipating the RG com- 
putation, we can arrange these perturbations into eigen- 
operators of the RG transformation as 
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(25) 



Using Eq. 



(4-d) 

we can compute the exact (3 function for 



dg 

~dl 



(4-d) 9 3 



(26) 



Also, from Eq. (|23|) . we obtain for the anomalous dimen- 
sion of the composite field $: 



(27) 



For d < 4, the RG flow equation Eq. (|26|) has an 
infrared stable fixed point at g* = \J(A — d). At first 
sight, this may seem puzzling because in Section [TT] we 
described the Feshbach resonance by an infrared unstable 
fixed point. However, this is just a matter of our restric- 
tion here to a particular critical manifold. We will see in 
Section llll B[ where we consider the full parameter space 
of the theory, that there is indeed an additional relevant 
operator here which corresponds to the relevant pertur- 
bation of Section HU and this is the detuning away from 



drd d x< 



v4v>2 



+ <s$ t $-/i(v4Vi -i>li>2) \ (29) 



We can interpret the three relevant perturbations as the 
chemical potential, /i, a measure of the detuning away 
from unitary scattering, 5, and the 'magnetic' field, h, 
which breaks the symmetry between ipi and ip2- It is a 
simple matter to relate S to the detuning v in Eq. (fTj). 
The T matrix for scattering between ipi and ■02 is equal 
to 5 2 G$. We have to evaluate this on-shell, for incoming 
particles with momenta k\ and fc 2 , which then requires 
G^(fc = k\ + k 2 ,uj = -i(k 2 /(2m 1 ) + fcf/(2m 2 ))) in the 



presence of a non-zero 5. Using Eq. (|22p in d = 3, and 
comparing the result with Eq. (TT]), we conclude 



27r(mi + m 2 ) 5 



ra\m 2 



(30) 



The RG flow equations for the /x, v, h perturbations are 
computed using the standard operator insertion methods. 
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From this we find the exact flow equations 
da 

i = 2 » 

dv ,„ ... 
^ = (2 - V(9)> 

n = 2h 



(31) 



Using the fixed point value of 77 in Eq. (f28[) . we see that 
the fixed point values of these scaling dimensions are 
identical to those of the three relevant perturbations ob- 
tained in Section [IT] for a field theory of atoms alone. 

Universal correlations of the finite density quantum liq- 
uid can now be computed as an expansion in the non- 
linear coupling g = g* . Because g* 2 is of order (4 — d), 
this leads to an expansion in powers of (4 — d). Indeed re- 
sults from such an expansion have already been presented 
in Ref. H 

As we noted in Section U the resulting (4 — d) expan- 
sion does have some limitations. All fluctuation correc- 
tions to the bare action S c + S p are of order (4 — d) , and 
are consequently unable to change the mean-field ground 
state obtained when [i, v, and h are of order unity. This 
has the consequence that a direct expansion in powers 
of (4 — d) is only able to describe the fully gapped su- 
perfluid, or the normal ground state with a single Fermi 
surface. States with 2 Fermi surfaces, or normal states 
at T > are only accessed after a redefinition of parame- 
ter scale s 25 ' 27 , and are more conveniently accessed by the 
method introduced in the following section. 



IV. l/N EXPANSION FOR THE 2N 
COMPONENT FERMI GAS 



Following, Ref. 130 . we consider fermions tp a with a = 
1 . . . 2N which have an Sp(2A) invariant interaction with 
the boson <E>. The critical theory for this case is (for 
simplicity we take the fermions to have equal mass m): 



(Jo 



8 



V 2 

Am 



roc $ 



&J af) i> a i>f3+H.c.)\ (32) 



where J a ^ is the Sp(2A) invariant tensor 



/ 1 

-1 



J 



\ 



(33) 



which is the generalization of the antisymmetric (4— d) a/3 
tensor for SU(2) = Sp(2). It is easy to work out the exact 
[3 function for this case, and we find 



dg (4-rf) 3 

de = ^ 9 - N9 



(34) 



Note that the fixed point value of g = g* is of order 
l/N. This implies that a l/N expansion is possible at 
any value of (4 — d) . 

Formally, the l/N expansion is obtained by dropping 
the bare quadratic terms in the critical theory for $. This 
is because, as we have seen in Eq, (f2"2"j) . the self energy 
contributions to the <E> Green's functions are more singu- 
lar than the bare terms, and dominate in the universal 
scaling limit. 

It is now possible to explicitly write down the theory 
which generates the universal properties to all orders in 
l/N. Henceforth, all momentum integrals are implicitly 
evaluated in dimensional regularization, and so we set 
roc = at the outset. We use the action S c + S p in its 
Sp(2A) invariant form, rescale $ — > &/go, and finally 
integrate the tj) a fermions out completely. This yields 
a theory for the single complex field $ alone, with the 
effective action 



iS<j> = 



\ 9oI ddxdT ^( x ' T ^ 



Trln 



or \ 2m 



+ Re[$]r 1 +Im[$]r 2 



(35) 



Here we have rescaled S — > N8, and r 1 ' 2,3 are the 2x2 
matrices in Nambu space. With the explicit factor of 
in front of the effective action for a l/N expansion 
can be generated for all physical properties. 

The above approach may not appear very different 
from the usual fluctuations-about-BCS theory of the 
paired Fermi gas. However, our formal key observation, 
which follows from the existence of the RG fixed point 
described in Section |H] and IIIH is that upon evaluation 
of all integrals under dimensional regularization, we au- 
tomatically obtain the scaling functions associated with 
crossovers across the Feshbach resonance. 

The effective action (|35p leads immediately to an ex- 
plicit N = 00 result for the grand potential which obeys 
the universal scaling form in Eq. In principle, this 
result can be obtained for an arbitrary spatial depen- 
dence of 4>, which would then allow for the spatially 
modulated 'FFLO' phases (magnetized superfluids with 
Cooper pairs at a finite momentum kp^ — kp^ determined 
by the Zeeman-split Fermi wave- vectors kpf and kp[)2^; 
this extension is discussed briefly in Appendix [A] We 
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restrict our attention here to a constant $, in which case 



TV 



47T 



1*1'- 



(2tt) 3 



Tln ^ + e -(\/(P 2 /(2™)-t>) 2 + m 2 -h)/T^ ^ 
+Tln (l + e -(V(P 2 /(2m)-M) 2 + l*l 2 + >>)/T^ 




i*r 



2m 



The superscript on J 7 denotes that expression above is 
divergent in the ultraviolet. As discussed above, the mo- 
mentum integral has to be evaluated in dimensional reg- 
ularization, by a process of analytic continuation from 
a convergent result for d < 2 to d > 2. We do this 
by adding and subtracting m\$\ 2 /(p 2 + A 2 ) from the 
integrand, where A is an arbitrary momentum scale. 
Then, the subtracting integrand is explicitly convergent 
for d = 3, while the compensating term can be evaluated 
analytically by dimensional regularization. This yields 



T _ (mv mAT(l - d/2) 
N ~ I 4tT _ 



1*1 



d 3 p 



( 47r )d/2 J i - 1 J ( 2 7t) 3 
T \ n (l + e -(^/(p 2 /(2rn)-^ + m>-h)/T^ (37) 
. T \ n fl + e -(V(P 2 /(2m)-M) 2 +l*l 2 +ft)/^ 




|$|2_ 



P 

2m 



m|$| 



A 2 



A ^ 



d d k 

(27T)° 



where the dimensionally regularized term can be evalu- 
ated directly in d = 3. It can now be checked that the 
above expression is independent of A, and is most easily 
evaluated at A = 0. The procedure above also shows that 
the subtraction is only needed for the |$| 2 term, but not 
for any of the higher powers in $. 

The phase diagram is determined by minimizing 
Eq. (|3"T|) with respect to <&. The resulting free energy 
obeys Eq. |J4]), and so yields a universal phase diagram as 
a function of fi, h, v and T. We will present results for as- 
pects of this phase diagram, along with 1/N corrections, 
in the subsections below. 



A. Superconductor to normal transition with 
increasing T 



At sufficiently high T, it is always the case that (*) = 
0. In this case, we can expand <S$ about <f> = for an 
arbitrary spatial and temporal dependence of To first 
order in 1/N, we need to expand the action for $ to 
fourth order: 



K 2 (k, w n )$t (fc, w„)*(fc, uj n ) 



^n( T E / ^^4(fc i! ^ n )$t(fe 1)a , ln )$t( fc2)W2n)$(fc 3 )W3n) 



X $(fci + k 2 - fc 3 , UJi n + UJ 2n - ^3n) 



Here 



K 2 {k,u n ) 



mv r d d p (1- f[{p + k) 2 /(2m)- fi-h]- f\p 2 /{2m}- fi + h] m 
4^ ~ J (2ir) d \ -iu n +p 2 /(2m) + (p + k) 2 /(2m) - 2^ p 2 



(38) 



(39) 



where /(e) = l/(e e / T + 1) is the Fermi function, and we have performed the dimensional regularization with the 
subtraction of the m/p 2 term. Also, we will only need K4 when two of its arguments have zero momenta and 
frequency, in which case it equals K4,{k,u> n ), given by 



Ki(k,u) n ) 



1 f d d p 

2 J (2n) d 



I ((k - p) 2 /{2m) -fi- h,p 2 /(2m) -p + h, 0J n ) 
+ I((k- p) 2 /(2m) -n+ h,p 2 /(2m) -fi-h, uj n ) 



(40) 
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where 

tanh(ei/(2T)) + tanh(e 2 /(2T)) 



7(ei,e 2 ,w n ) 



4e 2 (ei + e 2 - iui n ) 2 
tanh(ei/(2T)) + tanh(e 2 /(2!T)) - (e 2 /r)sech 2 (e 2 /(2T)) 

863 (ei + e 2 - «Wn) 
tanh(e 1 /(2T)) - tanh(e 2 /(2r)) 

8e^(ei - e 2 - iw„) 



Now we lower the temperature to critical temperature 
T = T c at which there is an onset of superconductivity. 
To order 1/iV, the condition for T c is given by 



tf 2 (0,0) 



2T 
~~N 



E 



(2ir) d K 2 (k,u; n ) 



(42) 



We will restrict our analysis of the consequences of 
Eq. (|42l) to the balanced case at unitarity, h = 0, v = 0, 
for which precise Monte Carlo results are also available. 
This will allow us to test the accuracy of the 1/N expan- 
sion. 

We tabulate the results of some integrals needed for the 
1/N computation. At N = 00, the transition is obtained 
by the condition 



were obtained explicitly. The summation over these 
leading terms was performed analytically, along with 
the required time-splitting convergence factors required 
for canonical bosons. After subtraction of these lead- 
ing terms, the summation over the remaining subleading 
terms was performed numerically, and converged rapidly. 

Using the first two results in Eq. (|4"4")) and the condition 
Eq. (|4"2")1 we obtain the critical chemical potential: 



= 1.50448 



T=T C 



2.785 

N 



+ 0{1/N 2 ). (45) 



K 2 (0,0) = 



T=T a 



= 1.504476695 (43) 



by 

P_ 

N 



The total density of fermions, p, is given to order 1/iV 



d 3 fc 
8^ 



All numerical values tabulated below are for the value of 
p in Eq. |g3J) and for h = 0, v = 0. 



(46) 

As is conventional, we relate p to a "Fermi energy", ep, 
by 



T E 



_ dK 2 (o,o) 

dp, 

d 3 k K 4 (k,ui n ) 
8tt 3 K 2 {k,uj n ) 

k 2 
if ' 
1 J 



d 3 fc 



8?r 3 V 2m 



M = 



2 ± /"Of - 

dp J 8tt 3/ \2m 

2T [ ^.ln(l+e-^ 2 / 2m -^/ T ) 
J 8^ 



0.018671(2m) 3/2 T~ 1/2 
0.0263(2m) 3/2 T 1/2 

0.096549(2mT) 3/2 

0.056179(2m) 3/2 T 1/2 

0.2258(2mT) 3/2 

0.13188(2m) 3/2 T 5/2 
0.1357(2m) 3/2 T 5/2 

(44) 



e F 



(3ir 2 p/N) 2 / 3 



2m 



(47) 



Using the third, fourth, and fifth equations in Eq. (|44j) . 
the value of p in Eq. |45|) and the relations Eq. (|46|) and 
(14"T)) we obtain 



T 



= 2.01424 



T=T C 



5.317 

N 



0{1/N 2 ) 



(48) 



Finally, the pressure P is the negative of the grand 
potential, and so at order 1/N 



— = T 

N 



d 3 fc 
8^ 



21n(l + e-( fe2 / 2m -^/ T ) 
T 



N 



^ln [Mk,u n )] 



(49) 



Using the third and the last two relations in Eq. (|4~4"1) , the 
value of p in Eq. (|4"5|) , and Eq. l|49p . we obtain 



We make a few remarks about the numerical tech- 
niques used to obtain the results in Eq. (|44|) . All fre- 
quency summations were evaluated on the imaginary fre- 
quencies. For all of the expressions, the leading terms at 
large uo n , such as terms of order ln(oj„), l/u> n , 1/w 3 / 2 , 



P/N 



(2m) 3 / 2 T 5 /2 



T=T C 



0.13188 + ^j-^+O(l/AT 2 ) (50) 



Although the 1/N corrections in Eqs. IjlSjl. (|3S]l and 
([Sn]), are quite large, direct evaluation of the expressions 
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at N = 1 yields values that are in reasonable agree- 
ment with recent Monte Carlo results 3 ^, which obtained 
fi/T c = 3.247, e F /T c = 6.579, and P/(2m 3 / 2 T c 5/2 ) = 
0.776. The leading N = oo contributions in the expan- 
sion are the contribution of fermion loops, while the 1/N 
corrections come from the boson loops: the distinct phys- 
ical origins of these contributions indicates that it is not 
meaningful to compare their values as a test of the accu- 
racy of the expansion. All subsequent 1/N contributions 
are corrections to either the boson or fermion loops, and 
the agreement with the Monte Carlo results is an encour- 
aging signal that such higher order corrections are indeed 
small. 



1. Relationship to other work 

The structure of our leading 1/N corrections bears 
some similarity to the recent work of Haussmann et a/ 37 . 
They consider the single fermion loop (our N — oo terms) 
and the single boson loop (our 1/N corrections) at an 
equal footing. Thus their approach can be viewed as 
effectively similar to a self-consistent 1/N computation. 
However, their approach breaks down near T c , and leads 
to unphysical results. Near T c , the value of ^(0,0) ap- 
proaches 0, and is also allowed to become negative in 
a self-consistent approach. Indeed, a direct attempt to 
obtain T c by solving Eq. (|42p self-consistently for arbi- 
trary N leads immediately to this difficulty: we obtain 
if 2 (0,0) < already for T > T c , and the fluctuation 
propagator in the second term has a negative 'mass' and 
so is obtained by performing a functional integral with 
an unphysical inverted Gaussian weight. However, orga- 
nizing the perturbation theory order- by-order in 1/N (as 
we have done) avoids this difficulty, and leads to a well- 
defined and controlled expansion for the thermodynamic 
properties at T c , with all boson loop propagators always 
obeying Re[K 2 (k,uJn)] > 0. 



B. Phase diagram at zero temperature 

We will now present some results on the universal 
phase diagram at T = 0. This is determined by 3 
parameters: /i, and h. However, these parameters 
are dimensionful, and the phase diagram can only de- 
pend only upon two independent ratios of these parame- 
ters which have zero scaling and engineering dimensions. 
From Eq. ([4]) , we can deduce that the phase diagram de- 
pends upon the ratios of fi, h, and v 2 /(2m). The sign 
of h is immaterial, while those of /i and v are certainly 
significant, and one or more of these parameters may be 
zero; consequently we cannot explore the complete phase 
diagram in a single two-dimensional plot. 

First, let us begin with some general considerations at 
h = 0. As long as the detuning, v, is positive, it will 
pay to have particles in the ground state of the grand 
canonical ensemble only for /i > 0. However, for v < 0, 



the negative detuning implies there is a molecular bound 
stated at energy —v 2 /m. So we can expect a non-zero 
density of particles for 2/i > —v 2 /m. This reasoning led 
to the phase diagram shown in Fig.O We will see shortly 
that this threshold value of \x is indeed obtained in the 
N = oo theory. 

For /i / 0, we will limit our analysis here to phases 
with ($) spatially uniform in the N = oo theory. At 
N — oo and T = 0, the free energy density in Eq. (137)) 
at A = reduces to: 



T 

N 




where h > was assumed without loss of generality. At 
h = 0, we can check for the transition in Fig. [5] by ex- 
panding Eq. (jo"Tj) to order |<I>| 2 for \i < 0. This yields 



T_ 

N 



^ = 1*1 



m,v 




_ 4tt j 


1 8tt 3 I 



2(p2/(2m) + |m|) P 



(52) 

Evaluating the integral, we observe that the co-efficient 
of |<i>| 2 changes sign at fi = —v 2 /(2m), thus establishing 
the location of the phase boundary in Fig. [2] at TV = oo. 
It is also easy to see that there are no corrections to 
the position of this phase boundary to all orders in 1/N: 
the <E> propagator has spectral weight only at positive 
frequencies at the phase boundary (including a pole at 
lo = k 2 /(4m)), and this implies that all loop corrections 
vanish. 

Let us now turn to h ^ 0. The superfluid phase in 
Fig. [2] has a gap, equal to |$| ('spin' excitations), and 
consequently a non-zero \h\ < |$| has no effect on the 
superfluid phase. The density of the particles remains 
balanced in this regime, and depends only upon the val- 
ues of n and v. 

We show a finite h phase diagram in Fig. [3j obtained 
by determining the global minimum of the N — oo free 
energy in Eq. ([51} with respect to variations in a space- 
independent <E> — a remarkable amount of information 
emerges from the minimization of this single function. 
The structure of this universal phase diagram has similar- 
ities to that obtained in the resonance- width expansion 14 . 
Because the physics is invariant under the change of sign 
of h, we plot the phases as a function of the dimension- 
less parameters n/\h\ and v / \Jlm\h\ In addition to the 
conventional gapped superfluid state with \h\ < |$| and 
balanced densities, there is also a more exotic superfluid 
phase in the phase diagram. The Luttinger theorems de- 
scribed in Ref. [38l . require that any ground state with 
an unbalanced density must have at least one Fermi sur- 
face (as long as the 'spin' rotation symmetry about the 
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FIG. 3: (color online) Universal phase diagram at T = 0, 
ft 7^ 0, and N = oo. The axes have been scaled by to 
the dimensionless parameters fi/\h\ and u/ y/2m\h\. The den- 
sity in the superfluid is balanced, except in the shaded (red) 
region representing a 'magnetized superfluid', which also has 
a single (1 xN) Fermi surface of (say) up spin fermions (for 
h > 0). The 1FS and 2FS phases are non-superfluid states 
with N and 2N Fermi surfaces respectively. The thick line 
is a first order quantum phase transition, while the thin lines 
are second order transitions. 



'field' h, pointing in the z direction in spin space, remains 
unbroken) . The shaded region in Fig [5] is such an unbal- 
anced superfluid phase, which has a single Fermi surface 
of one of the fermion species (that is, N rather than 2N) 
co-existing with the superfluid condensate. Modulated 
FFLO phases^ have (<&) ^ and space-dependent, and 
have Fermi surfaces: we briefly discuss such states in 
Appendix but have not yet undertaken the numerical 
analysis to determine their structure and stability — we 
will do so in future work. 

Other views of the same N = oo phase diagram are 
presented in Figs. [4] and [5] The only other stable states 
we find in our TV = oo theory (with $ space-independent) 
are the normal states with ($) = 0. This appears across a 
first-order transition indicated by the thick line in Fig. [4] 
As usual, the existence of this first-order transition as a 
function of chemical potential means that a system with 
fixed total density will undergo phase separation. In the 
present N = oo theory, one of the phase separated phases 
will be a superfluid with no population imbalance (fluctu- 
ations are not included at N — > oo). The other phase will 
be one of the normal states in Fig. [U which are present 
for h > h c (v). This normal state has either N or 2N 
Fermi surfaces, depending on whether the external field 
is strong enough to fully polarize the fermions. The criti- 
cal field at the Feshbach resonance is h c (0) — 0.807125^, 
and in the range h c (Q) < h < /i the normal state has 2N 
Fermi surfaces. 

It is also interesting to note that the local minimum of 
free energy at <& ^ survives at fields h as high as h c i(y) 
shown in Fig. 2] as the dashed line. This local minimum 
of free energy is global only for h < h c (v) < h C 2{v). 
Similarly, the free energy has a local minimum at <f> = 
for h > h c i{y). Therefore, metastable superfluid and 
normal phases can exist in the range h c i{v) < h < h c i(y). 




FIG. 4: The same universal, zero temperature phase diagram 
as in Fig. [3] but for \i > 0. Now we have scaled the axes 
by and they measure detuning and the field respectively. 
The first order phase transition between the superfluid and 
normal phases occurs at h — h c (v), plotted as the thick line; 
h c (0) = 0.807125/^. The dashed faint line denotes the "up- 
per critical field" hcziy) (equal to the fermion BCS pairing 
gap), below which the superfluid may be found at least as a 
metastable state. Similarly, the dotted faint line denotes the 
"lower critical field" h c i(y), above which the normal state 
may be found at least as a metastable state. 
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FIG. 5: (color online) As in Fig. [4] but for fx < 0. The shaded 
(red) region is a 'magnetized superfluid'. 



We now consider \/N corrections to the equation of 
state of the normal regions of the phase diagram. The 
pressure is calculated from Eq. (|49|) . which at zero tem- 
perature takes the form: 



P 

N 



/i(2m/i)5 
15tt 2 



1 - 



I' 



1 - 



h\ ( h 

1 + - 1 + - 



SP_ 

N 



(53) 



The calculation of the 1/N correction 6P is presented 
in the Appendix IBl These 1/N corrections are non-zero 
at T = only in the region where both fermion species 
are present in the ground state i.e. in the 2FS state. 
This state appears only for \x > 0, and so we assume a 
positive (i below. The scaling form in Eq. ^ implies that 
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FIG. 6: (color online) 1/N correction to pressure in the nor- 
mal state (h > h c i(v)) at zero temperature for fi > 0. In the 
limit N — > oo the faded pink portions of the curves should be 
discarded since they lie within the superfluid phase. 

at T — this 1 /N correction can be written in d = 3 in 
the form (/i > 0) 

where Fgp is a universal function. Our results for the 
function F$p are derived in Appendix [B] and plotted 
in Fig. [5] as a function of h/fi at various values of de- 
tuning from the Feshbach resonance. Virtual pairing 
fluctuations tend to increase pressure, and more so for 
longer lived pairs. However, for h > /i the fermion gas 
is fully polarized and no pairing fluctuations remain to 
contribute the increase of pressure. Note that the shown 
SP is valid only when the superfluidity is destroyed. This 
happens for h > h c (v), which in the case of infinite N 
corresponds to the dark portions of the curves in the 
Fig. El 

The results in Eq. ((53]) and Fig [6] can have a number 
of important consequences. Among the most important 
is that the position of the phase boundaries in Fig. Q] 
will be moved significantly, especially for negative detun- 
ing; additional phases not present at N — oo could also 
appear. However, determining this change requires de- 
termination of the 1/N correction to the energy of the 
superfluid state: this is a computation similar to that in 
Appendix [B] but more involved because of the broken 
symmetry in the superfluid state. This computation we 
also defer to future work, as we discuss further in Sec- 
tion |V1 Also, the equation of state of the normal state 
implied by Fig [6] has significant consequences for recent 
experiment a 28 ! 29 , and these will be described in a sepa- 
rate publication. 

V. CONCLUSIONS 

This paper has presented a unified understanding of 
quantum liquids near unitarity in various dimensions us- 



ing a renormalization group approach. The universal be- 
havior is associated with RG fixed points describing the 
interactions of a few 'atoms' and 'molecules'. This fixed 
point of the zero density theory is^ nevertheless a useful 
and powerful starting point for understanding properties 
of the quantum liquid with a finite density of particles. 
We view the RG fixed point as a description of a quantum 
critical point obtained as the chemical potential of the 
particles is varied. The quantum (multi-)critical point 
then describes a quantum phase transition (in the grand 
canonical ensemble) between a vacuum ground state with 
zero density of particles, and a finite density quantum 
liquid; see Fig [2] We have shown how this perspective 
allows us to bring the standard and well-developed tech- 
niques of critical phenomena to the computation of the 
phase diagram of such liquids. 

Our principal application of this method was to the 
attractive Fermi gas in d = 3. The fixed point theory 
can be described in either a (d — 2) or a (4 — d) or a 1/N 
expansion (for models with Sp(2A) symmetry). We ar- 
gued that the 1 /N expansion was best suited for the pur- 
poses of describing the onset of the non-superconducting 
ground state with increasing temperature or density im- 
balance, and presented results to order 1/N. We also 
obtained 1/iV contributions to the equation of state of 
the normal state. 

A number of additional results can be obtained in the 
1/N expansion, and these we hope to address in separate 
publications: 

• Single particle Green's functions can also be com- 
puted in the 1/N expansion by standard methods. 
In the normal state, these will yield information on 
the consequences of pairing fluctuations. 

• The existence and properties of the modulated 
FFLO phases are, in principle, also described by 
universal scaling forms such as Eq. Q. The form 
of the N = oo free energy of such phases is dis- 
cussed in Appendix [A] and a numerical analysis of 
the large parameter space of this result is needed. 

• Apart from the existence of FFLO and other 
phases, the phase boundaries in Fig HI will also be 
shifted by 1 /N corrections, especially in the regime 
of negative detuning. To compute these, we need 
the 1/N corrections to the free energy in the uni- 
form superfluid phase. This will have an expression 
which generalizes Eq. (|B10|) to the broken symme- 
try phase, and can be computed in a similar man- 
ner. 

• The equations of state obtained herein, and as in- 
dicated above, are expressed in terms of chemical 
potentials. These can be converted to equations 
of state as a function of particle densities by stan- 
dard thermodynamic methods. The local density 
approximation can then be used to convert the 
equation of state into particle density profiles in 
the presence of an external trapping potential. 
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APPENDIX A: 



TV = oo THEORY FOR FFLO 
STATES 



This appendix generalizes the N — oo result in Eq. (|37|) 
to states with a spatially dependent 

We expect that the optimum solution for $(2;) will 
have the symmetry of some Bravais lattice, and so write 



$(2;) = $ G e 



iGx 



(Al) 



G 



where G are the reciprocal lattice vectors. First, we 
need the energy eigenvalues of the Bogoliubov-de Gennes 
equations in the presence of such a We can expand 

the eigenmodes (u(x),v(x)) also in plane-wave Bloch 
eigenstates 

u(x) = e*P' x u G e iG x , (A2) 

G 

and similarly for v(x). Then the eigenvalue equation is 
•(p + G) 2 \ _ 

fl I Ug T 2_^i ^G-G'"G' = ep(*)UG 



2m 

:p + g) 2 

2m 



G' 



«G + ^ $ G-G' U G' = e p (*)l'G, 



where e P ( ( I ) ) are the Bloch eigenenergies. We now have 
to choose a truncation of the set of reciprocal lattice 
vectors, and then numerically diagonalize the equations 
above to obtain these eigenenergies. An infinite set of 
positive eigenenergies will be obtained for each p in 
the first Brillouin zone. In an extended zone scheme, 
these eigenenergies can be rearranged to obtain a single- 
valued and positive function of p, e p (<I>), where p ex- 
tends over all real values in (i-dimensional momentum 
space. This function can be chosen so that for <E> con- 
stant, e p ($) = ((p 2 /(2m) - v) 2 + \<S>\ 2 f/ 2 . 

Now, it is easy to see that the generalization of Eq. 



n 4tt 4r 1 J ( 27 3 



G 



Tin (l + e -(<-(*)-'>)/ r ) + Tin (l + e -^ +h V T ) 



1 2m 



(A3) 



where the integral is over all momenta. As in Eq. (|36p this 
expression suffers from an ultraviolet divergence, and we 
cure this by dimensional regularization. At large p, the 
values of e p ( ( i ) ) can be obtained by perturbation theory in 
<£>, and the divergent pieces only involve terms to second 
order in <i>. Adding and subtracting this divergent term 
to Eq. (|A3j) , and analytically performing the dimensional 
regularization as in Eq. (|37|) at A = 0, we obtain 



T vnv ^ -s . . 



d 3 P 
(2tt) 3 



Tin (l + e -( £ p(*)- fe )/' r ) + Tin (l + e -(M*)+M/^ 



.(•)-f^-^-5D«»l 



P 



G 



(A4) 



It now remains to numerically minimize Eq. (|A4[) over 
the set of values of $ G . 



APPENDIX B: EQUATION OF STATE AT T — 

Here we present some details regarding the 1/N cor- 
rections to the equation of state in the absence of con- 
densate at zero temperature. These corrections appear 
only in the normal states with two (2 x N) Fermi sur- 
faces, so that we assume /i > 0. The equation of state up 
to the order of 1/N is generally obtained from Eq. (|4"5|) 
generalized to finite external fields h. 

At T = and in the absence of condensate it is possi- 
ble to analytically calculate K^. Consider first the imag- 
inary part of Eq. (|39[) after analytic continuation to real 
frequencies iu> — > uj + i0 + : 



1 



lm{K 2 (k,co)} 



(2tt) 3 
(p + k^ 2 



1- e h 



p 

2m 



6[uj + 2fj,- 



P 
2m 



2m M 
2 (p + k 



(Bl) 



2m 



Note that in the rest of this appendix we will denote by 
K2(k,u) the dependence on real frequency u>, which is a 
slight change of notation from Eq. (f3"5|) . Upon the shift 
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Im{K 2 } 




useful, but they can be applied only to the well behaved 
portion K' 2 (k, oj). Fortunately, it is easy to recognize that 
the portion of K 2 {k^oj) whose imaginary part diverges 
when oj — > oo (the second term in (|B3[0 is just a square 
root function: 



- S j--m{2^ + 0J + i^) = 



-2 w n oj, 



mp 
i 6 

47T 



(p 2 ) 



(B6) 



FIG. 7: (color online) A typical plot of \m{K' 2 (k,u))} versus 
frequency (here fj, = 1, m = 0.5, h = 0.2, k — 1). There 
are up to five kink frequencies, listed in (|B8[l . including the 
end-points. lm{K2(k,io)} is non-zero only in a finite interval. 



of variables p — * p k/2, the delta function fixes the 
magnitude of p: 



k 2 

p = p(k,Lu) = \/m(2^ + oj) - — 



(B2) 



while the remaining integral over spatial direction of p is 
carried out easily: 



lm{K 2 [k,u)} = lm{K' 2 (k,oj)} - ^-6 (p 2 ) 



where 

Im{K' 2 (k,Lu)} 

x R 

and 

R{x) = 



= (p* 

m(oj + 2h) 
pk 



(B3) 



(B4) 



R 



m(u> — 2h) 
pk 





(B5) 



The real part of K 2 (k,oj) can be obtained from the re- 
quirement that K 2 be analytic in the upper complex half- 
plane of oj. Hence, the Kramers-Kronig relations will be 



We now determine the real part of K' 2 (k,oj) using the 
Kramers-Kronig relation: 



oo 

Re{K' 2 (k,oj)} = - [ du'^-Im{K 2 (k,oj')} . (B7) 

IT J OJ — OJ 

— oo 

A typical plot of lm{K 2 (k, oj)} is shown in Fig. \7\ De- 
pending on the values of [A, h and k there are up to five 
frequencies where Im{K 2 (k,oj)} has kinks: 



OjQ 



W 2 



— -2^ 

4m 



(B8) 



0J Q 



0J 



0J 3 = 0J + 



- - - y/2m(n + h) 
m V 2 



\j2m( y \i — h) 



1 ( k 

- - + V2m(/x - h) 
m V 2 



oj 



- - + y/2m(fi + h) 
m V 2 



We evaluate the integral (|B7|) piecewise between these 
kink frequencies. Knowing both the real and imagi- 
nary parts of K' 2 (k, oj) for real frequencies, together with 
(|B3[) and (|B6[) . allows reconstruction of the full function 
K 2 (k,oj), analytic in the upper complex half-plane of oj: 
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,, ,, , vcw m k 2 . . 
K 2 {k,uj) = — + -^--m{2^ + uj) 



(B9) 



772 



k 2 

2m(/i + h) — — 



2p(w\) + p(u) In 



p(oj) 



In 



\pN-pM 



2m(p-h) - — 



2p{u> 2 ) In 



p(w) 



In 



Vp(w 3 ) -pM 



+ p(cj 4 ) - p(o»i) - 



2A: 

p(cj 2 ) -g(g) 
p(a> 2 ) +p(w) 



0J4 — lui + (u) — 2/i) In 



+ p(cj 3 ) -p(cj 2 ) - 



2k 



u>3 — LU2 + (u) + 2h) In 



U4 — UJ 
u>i — UJ 



Ll>3 ~ UJ 
U)2 - W y 



Here, p(oj) is given by (|B2| . where dependence on fc is 
assumed and omitted for brevity. 

Now we can calculate the 1 /N correction to pressure in 
the normal state at zero temperature. For this purpose 
we substitute (|B9) into (135)) : 



SP 



1 

N 



d^k 
(2^)3 



— t 

2tt 



In [K a (Js,iu)] 



integral above reduces to: 



5P_ 

N 



i 

N 



d J fc 
(2^)3 

— oo 
oo 



duj /'K 2 (k 1 uj-iO + ) 
2^ n \K 2 (k,uj + iO+) 



(Bll) 



(BIO) = jji^ J dk J duj k 2 &rg[K 2 (k,uj + iO + )] 



The change of variables ioj — > w allows the frequency 
integral to be evaluated on the closed path in complex 
plane shown in Fig. [5] Avoiding the branch-cut of the 
logarithm, and exploiting the analyticity of K 2 (k, uj), the 



In the second line we have used the fact that K 2 (k,uj — 
i0 + ) = K 2 (k,uj + i0 + ). This integral has to be evalu- 
ated numerically. The plot of pressure correction at zero 
temperature is shown in Fig. [6l 
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